NTOUW-01-07 

Exploring Skewed Parton Distributions with Two-body Models on the 

Light Front: bimodality 

o 

O : B. C. Tiburzi and G. A. Miller 

■ Department of Physics 

University of Washington 
Box 351560 



00 



Seattle, WA 98195 - 1560 
February 1, 2008 



> 

00 

0^ ■ Abstract 



We explore skewed parton distributions for simple, model wave-functions in a truncated two-body Fock space. 
Consideration of non-Gaussian wave functions is our main emphasis, from which we observe the distributions are 
bimodal (there can be two distinct peaks) in the Fock-space diagonal overlap region [x > £). We demonstrate 
that this behavior arises from convoluting two light-front wave functions to form the skewed parton distribution. 
Furthermore, the factorization Ansatze, which have been used often, are shown not to hold. At high t, we can write 
the skewed parton distribution in terms of quark distribution amplitudes and perturbative wave functions thereby 
providing new tests. 

^ ■ 1 Introduction 



Recently there has been considerable interest in the connection between hard inclusive and exclusive reactions, which 
has been, in part, due to the unifying role of skewed parton distributions (SPD's) 0, 0, pjj. Aside from being the natural 
marriage of form factors and parton distributions, SPD's appear when one tries to calculate hadronic matrix elements 
such as in Compton scattering || |], || , and the electroproduction of mesons ||, [t]] . 

There is much effort underway related to the measurement of these functions Intuitively clear, but simple 
models || [| [Io|, [llj have been used to provide first calculations (and more sophisticated approaches have been pursued 
||l2| , |r3|). In this paper, we attempt to gain intuition about the structure of these distributions, using simple light-front 
wave functions. We find a new feature, bimodality (or double peaks), which is rather general and thus independent of 
sophisticated formalism. To this end, we work only with simple, two-body, model wave- functions for scalar particles 
and express the SPD as a convolution of these light-front wave-functions. This is, of course, a (very) specific instance of 
the general light-cone Fock-space expansion of the SPD's carried out in |Q. Of crucial importance is the employment 
of non-Gaussian wave- functions. As we shall show, the bimodal nature of SPD's is absent for Gaussian wave- functions. 
We remark that the bimodality of these distributions has been encountered before in jl^] where the authors study the 
SPD's in the chiral quark-soliton model of the nucleon. There it was shown that contributions from discrete levels and 
the Dirac continuum interfere, resulting in a roughly bimodal quark distribution (for x > £). We show below that there 
exists a rather general explanation for the bimodal behavior in terms of light-front wave-functions. 

We review the definition and kinematics of the SPD in section || and proceed to the two-body Fock space in section 
|3|. Once in section |], we present the wave- functions used and then provide a visual array of SPD's illustrating their 
functional behavior. Finally we conclude briefly in section ^| 
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2 Definitions and kinematics 

The skewed parton distribution is a non-diagonal matrix element of bilocal field operators. Various conventions, reference 
frames, variables, etc. exist for the description of such an object. Below we adopt the conventions of the off-forward 
parton distribution of Ji p"5[ . This choice is natural for a symmetric treatment of initial and final hadron states and 
consequently time-reversal invariance will be manifest. 

We shall consider the SPD for a toy scalar "meson" of mass M consisting of two scalar "quarks", each of mass m. 
Such a model restricts the SPD to only one kinematical regime (one that is diagonal in Fock-space). This restriction 
in modeling SPD's was encountered from the start j[o|, and we take up the issue of extending model SPD's to all 
kinematical regimes in [jl6| . 

Let us parameterize the kinematics as follows. The initial meson momentum is denoted by P p = (P _ , P + , P ), with 
p ± = (po ± p3) and the final 

meson momentum is denoted P'* 4 . The average momentum between initial and final 
states is thus P p = (P p + P ,Al )/2. Defining xP + to be the momentum conjugate to light-front distance y~, we have jl7| 



F(x,U)= 2(P |+ + p, +) / %-j y ~ x P+ ( P ' I ^ (" r / 2 ) d+ /2) - d + ^ (- y" /2) 0(y- /2) | P ) , (1) 

where the skewness £ is defined by £ = — A + /2P + with A as the momentum transfer suffered by the meson, namely 
A M = P' M — P M , with P' ± = — Pj_ = Aj_/2, and the invariant t = A 2 . As spelled out in time-reversal invariance 
forces F(x, £, t) — F(x, — £, t). It is easy to see why this must be the case in Eq. (|l|) owing to the reality of F and the 
Hcrmiticity of the current operator. Without any loss of generality, we take £ > below. Notice the form of the above 
equation is such that when integrated over all x, we recover the electromagnetic form factor (calculated from the + 
component of the current operator) . 

Having set forth all kinematical variables, we must then use Lorentz invariance to derive any remaining relations 
among them. Since our final meson is still intact, P 2 = P' 2 = M 2 and thus P 2 = M 2 — t/A. Exploiting A • P — 0, we 
find A - = £P 2 /P + , and also an expression for the transverse momentum transfer 

A ±2 = -4£ 2 M 2 - t(l - £ 2 ). (2) 

This implies a maximal skewness £ TO = 1/ \J\ — AM 2 Ji for a given t. 

In order to unearth the physics hidden in Eq. ([!]), we insert the mode expansion Q of the scalar field operators: 

with the quark creation and annihilation operators satisfying the relation 

[a(k+,k ± ),a t ( K +, K - L )] = {2ir) 3 d(k+ - K+)S 2 (k ± - k^. (4) 

Evaluating the plus-derivative of the field operators and inserting this along with the mode expansion Eq. (^) into 
the expression for the SPD, we find 



F{x^t) = \j %gj k ^ K e(k+)9(K+)S(xP+ - gjt£?l)Ck+ + K +) ( P' I a t( K +, ^) a(k+, k^) | P ) . (5) 
4 J V2k+(27r)' ;i V2 K + {2ny ^ 

3 Two-body, truncated Fock space 

We now proceed with representing the SPD as a convolution of light-front wave functions. In order to do so, we must 
write the meson states as sums of their light-front Fock-space components. The light-front Fock-space is the ideal arena 
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in which to decompose hadronic states [ fL9| |20[ ], since, apart from zero- modes, the perturbative vacuum is trivial. In 
this formalism, physical matrix elements are then expressed as a sum of Fock-space component overlaps (in general 
these overlaps are non-diagonal in particle number). The general form of the SPD written in terms of iV-body Fock 
space components was carried out in Ref. |14| and is a natural choice for these distributions since, e.g. the positivity 
constraints |2l[] are automatically satisfied ]1C| . Since one has little knowledge concerning viable TV-body light-front 
wavefuntions, this formalism is limited to the lower Fock-space components. In order to gain initial intuition about the 
SPD's, we shall limit ourselves to the two-body sector at the cost of neglecting one kinematical region. 

With a truly empty perturbative light-front vacuum a|0) = 0, we can build our quark states from the vacuum in the 
usual fashion 

|fc)^|< Z ;k+,k ± )=at(k+,k i )|0). (6) 

We shall keep the Fock-space as small as possible by truncating at quark pairs. In this approximation, our meson states 
appear as 

|P)=2P+(2tt) 3 f (f[ d P!l_ dp ^ )6 2 >+(P-p 1 -p 2 )*(p 1 ,p 2 )\p 1> p 2 ). (7) 

Several points must be clarified about the above expression. We are using shorthand notation for light-front momenta. 
Upper-case letters refer to meson momenta, while quark momenta are denoted by lower-case. The components (+, _L) 
are now implicit and we have notationally condensed the three-dimensional light-front delta function into <5 2,+ . 

The light-front wave function ^ appears in Eq. 0. The remarkable property of light-front wave functions is that they 
depend only on the relative momenta of their constituents, not on the hadron's momentum |p2[ . Our two-body wave 
function is only a function of the plus- momentum fraction x, and the relative transverse momentum p^, namely 



*(pi,p 2 ) = il>{x = Z 1 + ,p- L = {l-x) pi -xpi ) (8) 

Pl+P 2 

^{pi,p 2 ) is a symmetric function of p\ and p 2 . We also choose the normalization 

It is worthwhile to state the normalization of the state | P ). The use of Eqs. (f|](]) gives 

(P3,P4 | p 1 ,p 2 )H^f{^ + {pi-p^ + {p 2 -p i )- + {l^2}) (10) 

Thus we have 

(P'| P>=2(P+) 2 * 2 >+(P'-P) / d2 ' +P f^ P2 S^(P-p 1 -p 2 )\^( Pl ,p 2 )\ 2 (11) 

J PiP 2 

At this point, we must separate out the meson's momentum by defining total and relative momentum variables: V = 
Pi + p 2 with (+,-L) understood, x = p\ jV + and p = (1 — x)p\ —xp^ . The variable differentials are related 
by d 2 ' + pid 2,+ p 2 = ^-cP^Vdxdp 1 - . Carrying out the change of variables allows us to perform the integral over d 2 ' + V 
trivially since now the delta function 5 2,+ (P— V) sits in the integrand. The remaining integral is over relative coordinates 
and is merely the wave function's normalization 



P' | P)=P + 5 2 '+(P'-P) 1 4^ 

J x(l 

= 2P + (27r) 3 <5 2 '+(P' -P). 



dxdp- 1 



_Ls|2 



|^(^,P^)| 



(12) 
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To evaluate the SPD, we insert the two-body expansion (0) into equation (JsJ) and calculate the matrix element 

< P < i i p > - 4P+ p^f^£L=i±^^mi^ (13) 

\Jk + h + (P + — k+)(P'+ — k+) 

Using the delta function to perform the integral d 2,+ n yields 

F( X , t, t) = p+ p'+ / ^^^^Eipr^ + p+ - fc < fc ^ + A > p - fc ) 



2(27r) 3 (l-a;)(x 2 -C 2 ) J Vl-£' 1 - £ 2 / V 1 + ^ 

where we see the result of the theta functions present in Eq. (||) is to restrict the value of ir: 1 > x > £. In the other 
kincmatical regime, — £ < x < £, the SPD represents the amplitude to annihilate a quark pair from the initial meson 
state. Naturally the SPD in this regime depends upon the higher Fock components, which we neglect here. 

The above result is consistent with the time-reversal invariance of our starting expression Eq. (|l|) since it is clearly 
even in £. In the forward limit (t = 0), we recover the toy model quark distribution function 

^°-^/ 2 ( 2 Xa^) ^ kl)Ng(4 (15) 

Finally in the case of zero skewness we should arrive at the electromagnetic form-factor (after having integrated over x) 
@. We find 

/rlrrlV^- 
— -f- ^(x, k^*0r, k 1 - +(1 - z)A x ) ee F(t), (16) 
2{2TT) ti X(l — X) 

which is indeed the Drell-Yan formula [^3| . As pointed out in Ref . jl(| , the restriction to one kinematical regime (x > £) 
does not allow for the sum rule to be verified — Eq. (16) is as much as we can hope for now. We shall take up the issue 
of verifying the sum rule in Ref. Jl6[ . 

Note that the form of the SPD does not support various factorization j4nsatee [^4] . For example, in the limit £ — ► 0, 
F(x, £, t) ^ q(x)F(t), though this Ansatz does satisfy the sum rule. Such a factorized form only holds for Gaussian wave 
functions when £ = 0. Counting the number of field operators in equation fll we see there is no physical basis for such a 
factorized Ansatz when £ = 0. 



4 Model wave functions 

It remains only to choose two-body light-front wave functions to explore the SPD derived in Eq. (|lj). Before we do 
so, it is wise to review the existing literature on modeling SPD's. Radyushkin investigated them for a scalar toy model 
H deriving analytical results. They were investigated for a bag model of the nucleon in in the chiral quark-soliton 
model [^2| and in 1 + 1 dimensional QCD [[[3). More in line with our above analysis, SPD's were modeled using light- 
front Fock-space components in |l(| (using TV-body Gaussian wave functions for low values of N), for the pion using a 
Gaussian wave function in|Tl]], and for the QED electron in pq] . 

We will use three different wave functions to contrast the behavior of the SPD: a Gaussian (G), the Hulthen wave 
function (H) of Ref. p6|] , and the (weak-binding) Wick-Cutkosky wave function (WC) of Ref. ^7j. Each wave function 
contains a normalization constant iV (determined from the normalization condition used above in Eq. (^)) and is chosen 
to be symmetric about x = 1/2. 

Our Gaussian wave function is taken to be 

fc( .,^)_2^»M.-sfeS, (17) 

TO 
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Figure 1: Quark distribution functions F(x,0, 0) plotted verses x for the models under consideration. 



where f(x) = x 2 (l 



x) 2 is similar to a valence quark distribution (the difference from the quark distribution of Eq. 

M 



is due to the mass term in the exponential). We choose — = 1 and — = 1.8. 
The Hulthcn wave function has the form 



my Nx(l — x) 



with ^ , ^ chosen to be the phenomenological deuteron parameters 



4x(l - x)a 2 + (2x - l) 2 m? + k 

;al deuteron parameters. 
Lastly the Wick-Cutkosky wave function in the weak-binding limit is 



I5(^-^). 



(18) 



771' 



Nx 2 {l-xf 



{Ax{\ - x)k 2 + (2x - l) 2 m 2 + k-L 2 ) 2 1 + \2x - 1| ' 



(19) 



with k = T^ma, the meson mass M 



2m — jma . We choose a suitably weak coupling of a = 0.08. The form of the 
wave function is basically Coulombic with a multiplicative retardation term. 

In Figure [l] we contrast the quark distribution function for each of the toy models under consideration. The models 
clearly differ in how narrowly the quark distribution is peaked. The Wick-Cutkosky distribution appears sharply peaked 
at x = 1/2 which is characteristic of a weakly bound system. 

Next in Figure ||, we use our Gaussian wave function ipc in Eq. (|l4| ) to reproduce a typical Gaussian SPD found 
in the literature. We note the distribution maintains its characteristic Gaussian shape as a function of £. The peak of 
the distribution, however, appears to be suppressed as £ increases up to £ m . Since (as we will show) the special form 
of the Gaussian wave function obscures the general structure of the SPD, we proceed to consider the power-law wave 
functions. 

Let us plot the SPD using Eq. @ for the Wick-Cutkosky wave function ipwc- In Figure [| we show the £ = 0.07 
slice of the SPD at —t/m 2 = 0.1. The distribution is bimodal, with the two peaks displaced from x = 1/2. To understand 
this feature, we crank up the momentum transfer to —t/m 2 = 100. The larger momentum transfer allows us to probe 
a larger skewness (since £ m increases with — t). At this higher momentum transfer, let us look at another £-slice of the 
SPD as a function of x shown in Figure |J. 

Perhaps the bimodal nature of the distribution is not too surprising, after all the SPD is a convolution of two hadronic 
wave functions. The weak binding Wick-Cutkosky wave function is by nature highly peaked at x — 1/2. The integration 
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Figure 2: Skewed partem distribution for 2-body Gaussian model using the parameters: m/ji = l,M/n — 1.8 and 
evaluated for -t/M 2 = 1.23. 
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Figure 3: Bimodality encountered in the Wick-Cutkosky Model's skewed parton distribution. The SPD is plotted as a 
function of x for —t/m 2 = 0.1 and £ = 0.07. 
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Figure 4: Bimodality in the Wick-Cutkosky Model's skewed parton distribution at high momentum transfer. The SPD 
is plotted as a function of x for —t/m 2 = 100 and £ = 0.1. 



over the transverse momentum in Eq. ( |l4| ) does not affect the location of the peak. Due to the convoluted nature 
of the SPD, the x location of the wave function's peak becomes a function of £. Let us denote the location of ^>'s 
peak — tp( a(a;,£),k ± ) with the arguments in Eq. jl^ ) — as x p (£) and that of ip* by x p (£). Using Eq. (|l4|), we find 

^(0 = ^(1-0 
^(0 = ^(1 + 0- 

Thus for fixed t, as the skewness £ is increased, the peaks of ip and ip* trek in opposite directions. For £ = 0.1, the peak 
of ip should be located at x p — 0.45 and that of ip* at x p = 0.55. Looking back at Figure^ shows us this is indeed the 
case. This allows us to understand the peaks' behavior for the full range of £ shown in Figure |5|, where we see a ripple 
which heads to smaller x as £ increases. Once £ > 1/3 we can be assured the SPD's peak is due to ip*. But the figure 
opens up new questions. What can we say about the height of the peaks? 

Looking at all the SPD's displayed so far for the Wick-Cutkosky model, generally the peaks at lower x which stem 
from ip aren't quite the lofty summits as those from ip!* which are at higher x. This asymmetry is due to the prefactor 
in the SPD: ^J^^fJ^) which generally increases with increasing x for fixed values of £. Thus the peaks due to ip* are 
preferentially enhanced while those due to ip tend to be flattened as they trek to lower x. 

Taking a second glance at Figure ^, we see another feature of these distributions. The peak which we attribute to ip* 
is amplified with increasing £. To be able to make definite statements, we take the asymptotic limit and appropriately 
analyze the SPD. (We will find that this leads to generally true qualitative information for lower momentum transfer.) 
The Wick-Cutkosky wave function has power-law dependence, and so we appeal to the approach of Brodsky and Lepage 
||l9f of approximating integrals over kj_ using the knowledge that our wave functions are peaked for low values of k -1 -. 
Looking at Eq. ( |l4[ ) and taking A 1 - very large compared to kj^ , we realize the dominant contributions come from either 
k^ w or k^ « — Aj_ for which one of the wave functions is greatly suppressed. This gives: 

F(x,Z,t/m 2 » 1) « 2{2n)3 X { ^ 2 { ) x 2 _ e) (^i)V>*(z2, (1 - x)A^ + 0*(:r 2 )V(si,-(l " , (20) 

where <f>(x) = J dk. 1 - ip{x,\t^) and we've abbreviated the plus-momentum fraction arguments as x\ = f^|, £2 = fzf 
and their average x — (x\ + x%)/2. Notice that the above equation reproduces the correct location of the peaks as a 
function of x. The Wick-Cutkosky wave function, as well as the Hulthcn behave as ip ~ 1/t 2 for large \t\. In a realistic 
calculation, one would use wave functions calculated from perturbative Quantum Chromodynamics (pQCD). These 
wave functions have different asymptotic behavior ("0 p qcd ~ 1 A) than our model wave functions since the dynamics of 
pQCD stems from vector exchange (unlike the scalar exchange implicitly employed by our models). The pQCD analysis 
of the pion's SPD was carried out in E8] . 
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Figure 5: Skewed parton distribution for Wick-Cutkosky Model, with the model parameter a = 0.08, evaluated for 
—t/m 2 = 100. This asymptotic momentum transfer allows us to see and understand the general features of the 
distribution. Notice we plot the natural logarithm of the SPD. 



Using the form of Eq. (|2C|), we can deduce whether the wave functions are being stressed by increasing £. The wave 
functions will be stressed depending on whether larger transverse momentum flows through them as £ increases. Thus 
we are concerned as to whether (1 — xjA 1 - increases with £ (for fixed t and x). Taking the derivative with respect to 
£ shows that the change in transversal momentum flow depends upon the sign of \t\{\ — £ 2 ) — (1 + £ 2 )4M 2 . Hence, 
if \t\/AM 2 < 1, then less momentum will flow through the wave functions as £ increases. Otherwise, there is initially 
increasing momentum flow (starting from £ = 0) followed by decreasing momentum flow with £ p^]. In cither case, as 
£ approaches its maximum value (the maximal skewness £ m )) there is less transverse momentum flowing through the 
wave functions in Eq. ( po|) and consequently the height of the summit (centered at x p ) will increase. 

This fact is evidenced by Figure ||. Given the (natural) logarithmic scale, we see the SPD increases by a factor 
~ 10 near maximal skewness. To quantitatively understand this rise in the distribution, we'd be wise to consider the 
transverse momentum flowing through the wave functions. The ratio of the transverse momentum flow (in asymptopia) 
to the mass is given by <5(£) = (1 — x)\A \/M. For fixed x (say near the peak of ip* as we go to maximal skewness) 
and fixed i, we can plot this ratio as a function of £ which is done in Figure [| Indeed near maximal skewness, there 
is a drop (~ 10 x) in S and consequently we expect a rapid rise in the distribution since (rather suddenly) the wave 
functions aren't as stressed. 

The arguments presented above have been general for power-law wave functions and so aptly apply to the Hulthen 
model wave function as well. We can find the SPD's for the Hulthen model by inserting ipH into equation |l4|. At the 
moderate momentum transfer of 1 GeV 2 the bimodality is visible in the £ = 0.15 slice of the SPD, see Figure 0. The 
peaks are not quite so distinct since the Hulthen wave function isn't as narrowly peaked about x = 1/2 compared to the 
Wick-Cutkosky. Consequently distinguishability between contributions from ip and tp* diminishes. As we go to higher 
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Figure 6: Ratio of transverse momentum flow to M: <5(£) = (1 — .-r)|A J -|/M, plotted as a function of £ (taken here for 
t/M = —20 and x = 0.9). Note the drastic drop as we approach maximal skewness. 
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Figure 7: Hulthen SPD for \t\ = 1 GeV 2 and £ = 0.15 illustrating the bimodality. 



momentum transfer, the trends become as clear as before. Figure || plots the SPD at \t\ = 20 GeV 2 which matches up 
to the Wick-Cutkosky SPD (Figure |^) in every detail. At 20 GeV 2 , the maximal skewness is £ m s» 0.77. In the plot, 
we have shown the full range of £. The bimodality is still visible for £ < 1/3 but the scale is dominated by the rise 
toward maximal skewness. We ran into to this in the Wick-Cutkosky model (and was the reason we previously used the 
logarithmic scale in Figure |). Taking our focus away from the bimodal peaks, we see the distribution rises (~ 5x) as 
£ approaches maximal skewness even against the competing prefactor 1 — £ 2 . The power-law Hulthen wave function is 
clearly recovering from being stressed as we approach maximal skewness. 

Finally we must remark, the altered view of Figure || not only allows us to view the bimodality more clearly, it also 
hides a feature of in our model wave function! If we could look into the large x shadows in Figure g|, we would notice 
a striking asymmetry present even at £ = 0. The cross section shown in Figure ^ highlights this feature. This built-in 
asymmetry persists in the SPD for £ 7^ |3(J . Recall that the SPD at £ = gives the xm-x- integrated electromagnetic 
form factor. In light of our previous work with the Hulthen model form factor p6| , we readily interpret the asymmetry 
as due to factorization breaking present in this model. The dominant contribution to the electromagnetic form factor 
in asymptopia comes from a peak in the near end-point region (x < 1) and dominates the contribution from the peak 
at x = 1/2 by a factor of In \t\. As \t\ is increased the asymmetry shown in Figure ^ will develop into a strong peak in 
the near end-point region giving rise to the logarithmic modification to the form factor. 

Of course logarithmic, asymptotic modifications to the form factor are also present for the Wick-Cutkosky model 
|pi[ . One finds a near end-point asymmetry for the Wick-Cutkosky SPD as well, however, due to the value of model 
parameter chosen (a — 0.08), one must go to considerably high momentum transfer (\t\/m 2 ~ 10 5 ) to see this effect. 

Having visited the power-law, model wave functions, we must come full circle. What about the Gaussian wave 
function we initially considered? We must remark that the logic leading to Eq. (|20| ) depends crucially on the fact that 
the wave functions have power-law momentum dependence. Thus arguments about transversely stressed wave functions 
simply do not apply and we cannot gain intuition about the SPD's amplitude. Moreover (owing to the peculiar nature of 
Gaussians) taking the convolution of two Gaussian wave functions yields another Gaussian wave function (certainly with 
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Figure 8: Skewed parton distribution for the Hulthen model plotted for \t\ = 20 GeV 2 . The trends at this high 
momentum transfer are now identical to those in the Wick-Cutkosky model. 
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Figure 9: Un-x-integrated Hulthen form factor for \t\ — 20 GeV 2 , i.e. the SPD evaluated for £ — 0. The factorization 
violating peak has already started to develop (manifesting itself as an asymmetry). 



a different effective plus-momentum, but a Gaussian nonetheless). Thus bimodality (which was due to the individual 
nature of ip and ip*) must disappear. Moreover, at large momentum transfer, we do not see an amplification of the SPD 
near maximal skewness. On the contrary, the Gaussian SPD's are suppressed. Indeed the behavior of Gaussian wave 
functions is peculiar (not just due to unimodality). Luckily we have a wealth of intuition about power-law wave functions 
and thus our understanding of the structure of skewed parton distributions need not remain cloudy. The amplification 
(with £) of the SPD is a possible signature for power-law dependence. Physical observables (the deeply virtual Compton 
scattering cross section, for example) depend on weighted x- integrals of the SPD. Thus maximal skewness amplification 
could be canceled by compensating behavior in the other kinematical regimes. 

5 Concluding remarks 

Above we have derived the form of the skewed parton distribution in terms of two-body light-front wave functions in 
the region x > £. We then used model wave functions in order to explore the resulting distributions and gain intuition 
about their structure. The key results are: Eq. fll4[), that the use of power-law wave functions leads to skewed parton 
distributions with two peaks, that factorization of F(x, £ = 0, t) of into separate functions of x and t does not hold, and 
finally Eq.([20|). This last item gives F(x,£,t) for very large values of t in terms of quark distribution amplitudes and 
perturbative wave functions. Thus SPD's give a new way to investigate perturbative treatments. 
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